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Abstract 


tt 1h sUiwti tliat the dynamical halo model off era a natural explanation 
for the form of the variation of the cixsalc-ray path length with energy. 

The variation above 1 CeV/nucleon can be underatood aa due to the variation 
of the dlffualon coefficient, and hence the resident time, with energy. The 
flattening of the curve below 1 lleV/nucleon la seen to mark a tranaltlon 
to a convection dominated regime wtiere the diffusion coefficient Is no longer 
the determining parameter. A fit to the observations yields a tialo outflow 
velocity of 8 km sec An attempt Is made to determine the overall scale 
of the halo and the dlff«islon coefficient using recent ^^Be flux measure- 
ments but the data do not agree well enough to pin di>wn ti.ese variables 
to within less than four or five orders of magnitude. 


1. lilt roiliu't ion 


It luiM boon long known timt coNmU* r<iyM truvol throutth S~t cm 
of Intornti’lliir mat tor h«‘foro onoapintt fri»m t!u* italaxy. In tho laal low 
yoara, howovor. It haa booomo ovidont that tho anxnint ot nviltor 
travoraod, or gramnutto, la not oonatant but la onor^y ilopoiulont . Abovo 
I lIoV por nuoloon tho moan loak;i)to path appoara to ilooroaao with «*nornv 
aa K whovo ci (.hilluaaon ot al. H72; Smith ot al I"*?!), 

Holovr I th'V por nuol«*t>n tho ovUli'noo la aomowh.it unoloar but If i*no 
foousoH onoa attontli«n on tho abuiuLinoo rat loa ot iIh'ho niioloi that aio 
lH*lii*\’oil to bo puroly sooomlary to tho paront nuolol thoir la ovldonoo 
tlwit tho moan loaka^o path la .1 oonatant ‘>-t* i;m om ^ . Tlila ovldonoo 
la not oonolualvo, howovor, ami varloua .lulhora havo dll forint; i*plnl»nia 
rogardliit; thla nuittor (Oimi*a and Krolor '.‘*78; I.o/nlak and Wobbor l‘*78; 
Kontoa it al . .‘)77). For tho proaont wo ahal I ailopt tho vlow tliat tho 
liuli’pinulanoo (or woak dopond.inoo) ot t tio coaralo-ray t;ranin.it;o bo low .1 low 
lIoV la roo I . Tho tollowliit; dlauoaalon mual thorofoio bo ro^ardi'd .la 
tontatlvo and aubji’ct to furthor Invoat li;at li»n ol tho v'baorv.it tonal 
a it uat Ion . . 

Thi’ boluivloi ot tho imMii loakat;o p.ith abovo I iJoV por nuoloon la 
us'ially dlaouaaod witliin tho I ratnowi'vk of tho "loaky box" mi’ilol with an 
o.soapo probability tli.it Inoroasoa with onorgy. Slnoo moat thoorloa ot 
chart;od partlolo propat;at Ion la dla rdorod m.i);not Ic tloKla prodlot .1 
diffusion ooofflolont that liuMoaaoa with partlolo onoitty (at riot ly 
aiioaklnt;, rluldltv) .in onorgy dopondont oao.ipo proKib 1 1 It v , or oqulva- 
lontly, roaldonoo t lim* la to bo oxpootod. 

For oxamplo, In tho oaao ol almjilo, ono dimonaional diffiialon In a 

» 

fi'nlon of 1 tno.ir alzo P tho moan roaldonoi' l inu* ti'r p.irtloloa ta t <il* /tC 
whor o < I a t h i* i! t f f ua i on o 00 ft lot ini t . 
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1110 nviln difficulty with Much a miidol 1 h Hn liuiblllty to explain 
Iho rolatlvo indopondonoo of the m«*an path lony,th on energy for partleloH 
below ~ I GeV/nucleon. There 1m no reason to expect the diffusion 
coefficltnU to exhibit this behavior therefore soim* additlomil ingredient 
is required. In this note I shall donxmatrate that the concept of a 
dynamical, outward flowing galactic h.ilo can supply this missing 
Itjgred lent . 

The notion that a tenous luilo about our galaxy could have a 
profund effect on the observed properties of cosmic rays is certainly 
not a new one, (See, for example, the refiTcnces in Stecker ami Jones 
1977 ) I'ntll recently, however, the halo has been regarded as a static 
diffusion volume that acted only as a passive storage region where cosmic 
rays could spend a great pt'rtlon of their lifetlnx's without encountering 
very much nutter It Is now appreciiited that it the nwiterlal of the halo 
Is moving it can carry the cosmic rays with it and this connective iik>t ion 
can s Ignit leant ly attect observed jvirameters of the cosmic rays (Joklpii 
197h, iVnes and Joklpii ld77n, b) , As we shall now see one i»f the parameters 
that can be affected is the mean path length traversed by the cosmic rays. 

In the ensuing dlscusston we shall consider a highly simplified 
nw>del of the galaxy. The disk of the galaxy that contains the cosmic 
ray sources and all of the nutter tluit contributes to thler grainruge 
will be represented by an inf lnlteslm.il ly thin disk of Infinite extent. 

Tlie halo will eUend In the X direction, perpendicular to the disk, to 
.an outer free iscape boundary at X • Tire mitflow in the luilo will 

be repri'sented by a cosralc-ray convection velocity, V in the positive 
X direction for positive X an»l In neg.atlve X direction for neg.it ive X, 
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Thiii 1m Mi'hinRat li'dl Iv 1 1 luHt ml tM it t lituro 1. Hu* aunli'l 

Im higitly It Imi'I it Its! in that it iKitort'ii aiiv variation of conmio-ray 
proitiict ion, or i^alaotlo paramotoiH ovor the diak of the );<ilaxy. Titere 
in evidence, however, tliat eoHmie raya do not travel about the diitk for 
dilttaneeN mtieh greater than I l,>e during their lifetime (Stecker and 
Jonea H??, Ornwa and Kreler l*>78^ therefore v»ne vanild exj vt that a 
local, one-dimenalonal nx«del wtnild be a reaaonable lepreaentat ion ol 
reality. rite diare)tard of the diak componenta flei»: thickneaa aa well 
aa poaaible varittiona of the convect Ive velocity with X are im*re aeriotia 
omiaalona. It ia expected, nevertheleaa, that reaulta obtained with thia 
initdel will ;iol ditter tiH* Kerlinialv from thoae obtained from nx're 
d»*i ailed mode I a 

In exa>n*.nln)t thia patticular nxulel Joklpil (l‘>7n) ahitwed that the 
mean a^e of the coNmic raya that are obaerved in the diak dependa ai);ntfi- 
cantlv on the aife ot the dlmenaionleaa panuneter \ ” M'/k. Tlila parameter 
beinjt eaaentiallv the ratio ot the convect ive velocity to the diffnaion 
ve’ocltv, la a ro<*aauie ot the importance v»f the outward convect ion, 

Ji'ktpll ah.«wed that the mean ajte or realdence time ot the comalc-ray 
part idea obaerv«>il in tlu' diak approached the limltln); lorma: 

T o'/iK, \ 'wv I tn 

T K/V’ . \ » ’ . 

Kiinat ion il> ia what one w\«nld expect lor a aimple diffnaion picture 
when convect Ion is unimportant. K^juat ion (2> ahowa, hviwever, that when 
conv«'Ct ion dominatea the mean a>;e beci'mea Independent ot the aire ot the 


CvMitalnment w'lnme 
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Thv presont «.uthur (Jonen 1978) ^ovc a simple explanation of this result 

* 

In terms of the "effective size" of the containment volume, D • Tlie 

concept of the effective size is readily understood If one considers 

a particle that finds Itself a distance x out in the expanding halo. 

2 

It will take this particle a time of order r ■ x /k to diffuse buck to 

the disk. During this time convection will carry it a distance outward 
2 

equal to Vt ■ Vx /k'. If this Is greater than the original dlst.mce 
It had to diffuse to return it Is very unlikely to mke It back. In 
other words if 


or 



(3) 


X > - tt n* (4) 

the particle Is effectively lost from the containment volurai’ Insofar 
as observations in the disk are concerned. 

Furthermore, using the concept of retrodlctlve probability (Jones 
1978) it is possible to show tluit those particles that are observed In 

the disk have spent their entire past history confined to a volume 

* * 

characterized by D rather than D. Of course. In the case D >0 (\nJ) 

* 

the actual boundary of the volume' prevails lutd D lias nei particular 

* 

significance. It Is straightforward to verify that substitution of D 

for D In ecpiatlon (1) yields equation (2) (to within a constant factor) 

* 

shwlng that 1) Is indeed the physically relevant confinement scale 
whenever It Is smaller tlun D, 

Turning our attention to the nx>an path length traversetl by the 
observed cosmic rays It Is evident that 


I “ VpT 


( 5 ) 


5- 


wh#r« < thtf nuMn p4th p tho ia*.in .tvitccr Jonslcy !*»>en bv 

Ch« cosmic rays, r Iho mr^n A|(e of cho cosmic r.tvs .tnd v Is tho spood of ch« 
particles. If the matter in the disk is of dcuMity p .ind confined to 
4 disk of half thickness «t the cosmic rays, confined to .1 voUime of h«lf 
thickness will see .t mean matter density 


p • • ih' 

TIm mean lifetime is itiven bv eipiation (1) iisinit 0^ instead of D 

* 

where P will be either P or P whichever is appropriate, Combining 
equations (5^, and (l> we obtain 

7 • -T vp a P' /\- . f7> 

y o 

We see now that if \ is large enough such that \ I and 
P* > P, P^ should be replaced bv P and the mean path length varies 
Inversely as <, If < is sixill enough for convection to domlaite, 
however, \»l and D' should be replaced bv P* ■ V and the mean path 
'ength is given hv 


f - 


_l 

\ 


vp a,V 




and is Independent of <, If the particles are still rel it 1 .’1st Ic at 
these energies f is independent of enorgv since v « c *’ur or non 
relativistic particles an energy dependence will entet rough the 
particle velocity v. 

“ • rhe P i f lus lo n-Convec t Io n Fquat ion 
rhe foregi'lng arguments are not precise and furthermore thev 
neglect the effect of the expanding !vilo on the particles energv 
Since we are interested in the energv dependence of the results It is 
Important not to neglev't energv in the derivation A ••.vre correct 



procedure Is to solve the equation for the cosmic-rav particle density 
NU.x.(.T) 


1 ^ vp 3N ^ V 3N :*N ,, , 

^ • 3 H vr ^ Tt ^ ^ 


o 


I**) 


where ■( is the path letmth, T is the particle kinetic energy and 
S is the source term that describes the particle injection. Fkiuation (*>) 
differs from the one employed by Jokipli (l'i7h> only in the addition of 
the tern describinti the chanite in path length. 

Wliat we desire is a steady state solution of equation (•<), A 
completely gener.al solution va'uld be quite difficult to obtain because 
the coefficients of the path length and diffusion terms are energy 
dependent hence the eqitation is not separable. However, if we assume 
that the pri.-nar>' energy dependence of the solution is due to a strongly 
varving injection term, eg. an inverse power law, we awiy ignore the 
energv dependence of these coefficients to obtain a reason.'.ble approxlnuite 
solut ion. 

If we adopt the following 


p(x) • Jp a6vx> 

S(\.T.O • 

o 

we obtain the so Ivit ion to eqvjatlon with 3N/3t ■ C> 


N 


:vp^a 


exp i- h' O 


a-exr[VU-D V\l> --V 
i 1-exp t-V:>'Tn 


where the mean path length is 


( 


^•co»'«[ l-expi-VT 1 

Vil -t- f l-exp(-iT> \->l} 


(10) 

(ID 

(ID 


(13) 
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If wi* Adopt for ttu' dlffuHlon coef f IcitMit K tho form 

1/2 

K - (15) 

whore ■ v/c and R iti tho particle rigidity, wo obtain an oxproHaion 
for 4 AH a function of rigidity; 

B 1/2 

L ;>ll-oxp(- 

< - , • (16) 

I + (T!-)( l-oxp(- ^ 

’ JJ r1/2 

In figure 2 wo Hhi*w a compilation of diita after OrmoH and 
Freler (1178) and a fit of exproHAioiiH (Ih) whore wo luive chosen the 
paramoterH 

_2 

L ■ p ac/V • 20gm cm 
o 'o 

I 1/'’ 

R Vn/K - ,'7 (I'V) 

o o 

We have also shown the effect of including a retarding potential 
♦ •' 220 MeV/nuc at I A.O. to approxinvite the effects of solar nk>dulatlon. 

Tlu* value of ♦ • 220 MeV/nucleon was found by Oarc ia-Munos et al. (H77a) 
to best order their iLita on cosmic-ray composition. 

Clioosing 2p a to be 5, lb M pc *■ In the local galactic region, 

U • 10 kpc , (ik>rdon and Burton, ll7b) we obtain V • 8 km sec ^ for the 
halo outflow velocity, TlUs velocity, while not large, has observable 
effects. For if tlie present Interpretat ion Is correct it is responsible 
for the turn over at low rigidities of the curve of cosmic-ray path 
length versus rigidity. 

3. Tlie Me.m Cv^smtc- Ray A^e 

We nkiy also determine that the ratio 


D 

K 


o 


X 10 (i.V) sec-cm 
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but to ilotormlno 0 or Boparatolv wo roqulro an adilttiomil ploco of 
oxpor Iraontnl data. Hu* mi*an a|(o of tho oormii' raya would bo auch an 
lndopondt*nt quantity and aworal oxpor Imonlora (Garc la-Hunos ot al l'»77b, 
Wobbor ot al l**77, IUt;on ot al 1*177, ButflU)tton *t al ld78) tuivo roportod 
valuoM lor thla qiwmttty l>a»od on nK'aHurt*mont a of tho flux of radloaotlvo 
^'^Ro la tho ooaralc-rav boain. I'no raiiat not think, lu'wovor, tiuit tho 
coamlo-ray a^oa roportod by thoso authora n*»y bo appllod dlrootly to 
our proaont propagation .u’dol. Aa w:ia jH'intod out by Prlahohop and I'tuakin 
(l''7M and Giitzburg and Ptuakln (1*^78) thoao agoa woro dorlvod from tho 
obaorvat Iona within tho hop^.y 'nooua or loaky box nx'dol of tho g;ilaxy anJ 
tho aamo obaorvallona oan load to qutto difloront agoa whon intorprotod by 
a diffualia or dil tualon-oonvoot ion nx'dol. 

To obtain mt'.uiinglul quant it ioa to apply to our iulhIoI wo must 
roturn to a px'ro priraltlvo oonoopt with roapoot to tho ^^^Wo nH‘aauronH*nta. 
Tho concopt ol surviving Iraotlon f la tv*t oonplotoly tu’dol indopondonl 
alnoo it roqulroa a oaloulatlon of tho flux of ^^^Bo that wv*uld bo 
obaorvod if it woro atablo against radliMctlvo docay. Aa wo aht>w in tho 
Appondlx tho path longth distribution funotion is what is noodod to 
comparo ono nx'dol to anothor i horntoro, if two dlfforont nx'dola 

prodlot tho aamo path longth distribution for atablo laotopoa thoy will 
prodlot tho aarao ratios of thoao atablo laot ,.a. Slnco tho loaky K*x 
nx'dol and tho proaont nxulol b**fh prodlot an oxponont lal path longth 
dlairlbution wo nuiy tako tho surviving fraotlon f to bo a aufflclontly 
primitive not ion that it nw»v bo usod without nxuUfication in olthor nx’dol. 

Although tho varlouN oxpor Imontors tiso tho aamo galactic propagation 
nx'dol thoy do not uso oxaolly tho samo calculation procedures, spallation 
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croaa ■•ctiona, or solar nodulacion nodala to compute the non-decay flux 

of and hence the surviving fraction* Therefore, solely for the 

sake of bringing the various experimental results into a common framework, 

we shall compare the reported values of the Be/ Be flux ratios to these 

calculated by Raisbeck and Yiou (1977) for the case of no decay of the ^^Be* 

10 9 

These authors calculate a value of 0*6 for the no decay Be/ Be flux ratio 
at 1 A.U* from the sun, virtually independent of energy* 

The surviving fractions are then simply given by 
f ■ (^^Be/^Be) obs./0.6. 

For propagation models that yield an exponential path length distribution, 
as does uur model, we have 


L/ ^ + o 

where is the mean path length neglecting decay of the particles, 

•(q is the mean path length including the effects of decay, and o is the 
collision probability per unit path length. 

Equation (9) may be easily modified to Include the effects 
of radioactive decay to obtain 
_ 2^Po» 1 2(v^l) 

V Ip 3 ^ 



(17) 


(18) 


where 


F 

X 


taah(^ ) 

6+tanh 



(1 + 


4k 

^ 

V^YTo 


and yTq 


is the relativisticly dilated decay time of the 



(19) 


2.16 Myr. 


V ^ 


I 
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Slncc* wtf have previouiily chontMi valucii for x, d.ca/V, and V onft 

oan inaort oquatlon (18) In ecjiuillon (17) and iiolv*> tho ronultlnr. 

oxprofltilon for 6 ^nd ht'nt'o K In tahU* 1 wo ithow tho varlouM i«'dol 

10 4 

par;imt'tor> that can bo deduct'd from tho reported Bo/ Bo ration. It 
In Imodlatoly obvlouM that the comaic-ray lifetlmoa are conaldorably 
longer than thoao deduced for tho leaky bt»x DK^del, In some Inatancea 
as much as tuv* ordora of magnitude longer. It la equ^illy ohvloua that 
the preaent experimental aituatlon doea not define these parameters very 
alkirply, the allt>wed values. Including errors, of some qiuntltlea range 
over four or nx>re decades. It appears that we must await a considerable 
sharpening of the experimental situation lu'fore we can have believable 
• 1 .18 for these Irap^^rtant paramc'ters At present the result of Stocker 
and Jones (1977) limiting the value of D*si 3 kpc Imposes as firm a 
constraint as any derived from observations of the flux of ^^Be. 

4 . Cone luslon 

We have seen tluat the variation of cosmic-ray path length 
with energy nwy bo understood within the ontext of a simple diffusion- 
convection nxidel with a dynamical galactic lulo. A fit to the observations 
of the various nx'del p.irameters gives rather directly an outflow velocity 
of ~ 8 km sec Tills velocity, while quite small. Is, nevertheless, 
responsible for the turn over at K>w energy of the cosmic-ray path 
length versus energy curve This fit also gives a value of D/x ■ 

1.9x10 ^sec-cm ^ for a rig Id I tv R ■ 2G\’ where we assume X ■ K 

o 




TABLE 1 




Galactic parameters 

derived from Cosmic-ray 

^^Be Measurements 



K.E. 

f* 

D* T 

K 

Authors’^ 

(MeV/Nucleon) 

Surviving fraction 


,,-27 2 
(10 cm -sec 


80 

.14 +.07 

125'^ 

-7.5 

+1.9x10" 
^ -334 

31 2*^^° 
-18.6 

145-200 

33 

-.20 

, --H6.8 
-1.6 

79 9"^^® 
-56.3 

5 

• -3.9 

200-245 

13 

• -.11 


-464 

48 9'*‘2*2x10 
-35.4 

250 

.62 +.27 


-9.6 

984'^-3 
• -.763 

100-1500 

.37 +.18 


61«« 

6.98««:? 


/ A-Garcla-MuiK»z et al. (1977b), a. Webber et al. (1977), C. Hagen et al. (1977), D. Buffington et al. (1978) 
10 9 

* The Be- Be ratio reported by these authors may be obtained by multiplying the surviving fraction f 
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KUIURK fAl’TIONS 

F'lKuro 1. 5rhrnuit li'ii ot ah«tU*l of ^Mlart tr ooHmlr r<iv prop<iK<it (i<n , 

Fluiiro 2. MtNisui «*nta*nt M ot path IcnKtt) va. rlKiiHlv t rora varloiin autliora 
t'ktmplloil hy I’rini'H aiul Fr«*t«*r with oiirv«* of I'xpioMalon 

(Ih). Soo toxt lor oxplan.it ton ot paramot oi'm . Oat.i aro 
Iromi op«Mi oli U'OM_, ot al. (l*)77>; aol l»l I'tioloa, 

('alttwoH (l‘>77); Caro ta-Munoc ot al. 0*>77a); »*pon 

li laiii^lo. Smith ot al . (l‘»7l); Sipian*, .Ini luaai>n (l‘*74>; 
oloaotl trtan^lo, Wohhoi , l.ofniak, .nut Klnh 11*171); liox<i^i«n, 
SIwipiro ot al . (l‘>7l>. 
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Appendix 

CoiiHlder a point In the matter diak of the galaxy. The denaity of 
part idea of the apeelea that Itave traveraed between { and f-Kl< gm/cm^ 
of mat: er, neglect Ing collialona, la given by the path length diatrlbutlon 
tuner Um Tt(f)d{. If ttu* particlea are injected at '( ■ 0 with a rate-denalty 
we may write down the value of l'|(U) provided th** aourcea are all in 
the matter diak. 

Clearly 

Pj(0) - Si dr (Al) 

where dt la tlw time required for a particle to traverae d< of matter, 

dt - (df/dt>“^ df - (vp^,)“^df 

thua Pi(0> - Sl(vp^,) It, In addition, the path length distribution 

ftmctlon ia an exponential function aa in the present nx'del or In the leaky 
U»x model then we must have 

I’lCfl - — ^ exp(--f/fi> (A2) 

''I'o 

where all of the dependence on the diffusion coefficient, distance to the 
boundary, decay time tor radioactive species etc. must be carried by tlie 
mean path length . 

It we now consider the ef tecta of collisions we obse* •• tluit if 
is the colllsliut probability per unit path length the probability tluit a 
particle which luis traveled a length f has not undergone a colllalini ia 
given bv exp(-Ojf). Since a collision will transmute a particle of species I 
into one or nx're particles of a different species the actual path length 
distribution for species I particles is given by 


1*1 (f) (actual) 


exp(-xjjf) Pj(f) 


(A3) 



miii the totMl dtMiMitv of MUi'h pdrtli'leM In the iliitW In Jimt 


o 

S, — -I 

• 0/^4 ^ ^1^ lA"*^ 

v»\, 

tf 1 m «tn exponent la I tunet li«n. 


It one i'onMlJei'H the expeotevl itenaitv ot rad ioaet Ive |«artlelea ineUidln); 


the etteetH ot Ueeav Np and Ignorlni; deeav the dltterenee will be entirelv 


dne to t Ita* dltferenee In < and henee the aurvivlnt; traetton 
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